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Experimental determination of an unknown quantum state usually requires several incompatible 
measurements. However, it is also possible to determine the full quantum state from a single, re- 
peated measurement. For this purpose, the quantum system whose state is to be determined is first 
coupled to a second quantum system (the "assistant" ) in such a way that part of the information in 
the quantum state is transferred to the assistant. The actual measurement is then performed on the 
enlarged system including the original system and the assistant. We discuss in detail the require- 
ments of this procedure and experimentally implement it on a simple quantum system consisting of 
nuclear spins. 

PACS numbers: 03.65.Wj, 03.67.-a, 05.30.-d 



I. INTRODUCTION 

Given the state of a quantum system, one is able to 
calculate the results of any measurement performed on 
that system. However, to determine the state from the 
results of measurements, one usually has to perform dif- 
ferent measurements that are not mutually compatible, 
using non-commuting observables. This issue is some- 
times referred to as the "Pauli problem" , since Pauli dis- 
cussed it in 1933 Since then, interest in this issue has 
continued, as it touches the fundamentals of quantum 
mechanics f2| More recently, it was also found to be of 
practical importance in quantum communication [3J IHIS] 
, e.g., in quantum cryptography and quantum key distri- 
bution [SJ[7]. 

If we consider an ensemble S of A^- level systems, its 
quantum state is described by a density matrix p in A^- 
dimensional Hilbert space, which requires N"^ — 1 real 
parameters for its complete specification. These parame- 
ters can be determined experimentally from the outcomes 
of a series of different measurements on identically pre- 
pared ensembles. Since quantum mechanical measure- 
ments with an observable U, with the spectral decompo- 
sition Cl — ^aPa generate at most A'^ — 1 indepen- 
dent probabilities, at least {N"^ -\) / {N -\) = N+l mea- 
surements with noncommuting observables are required 
to fully determine the unknown state p. 

Techniques for the reconstruction of the complete 
quantum state from a series of measurements are com- 
monly referred to as "quantum state tomography" [8l |9j 
ITU] . Different versions of such techniques have been pro- 
posed, with the goal of obtaining the best possible infor- 
mation about the unknown state while using the smallest 
possible number of measurements. Since the number of 
measurements must be at least TV -t- 1, the task is thus 
to determine an optimal set of + 1 observables (see. 
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e.g.,[TT]). A solution to this problem was given by Woot- 
ters and Fields They found that the observables 

should be chosen such that their basis states are evenly 
distributed through Hilbert space "mutually unbiased" ) . 
This choice assures that the data redundancy is mini- 
mized and the information content is maximized. 

As the simplest example, we consider the state of a 
spin 1/2. Its density operator can be written in the form 
p — ^{1 + s ■ a), where tr are the Pauli operators and s 
is a dimensionless vector of length < 1 that specifies the 
position of the state in the Bloch sphere. 

The simplest approach to determine this state consists 
of measuring the spin components along the x, y, and z- 
axes, yielding 6 possible measurement outcomes (see Fig. 
[l] a) . A minimal set of measurements only requires four 
such probabilities. They may be chosen as the proba- 
bilities for measuring the spin operator components in 
four directions that are oriented like the face normals of 
a tetrahedron [13] . 

While these approaches all require a combination of 
measurements with incompatible observables, it is also 
possible to obtain the complete state information from a 
single measurement performed on a larger Hilbert space, 
provided the state information is first redistributed into 
this extended space [H [TSl US] [T^ . Ref. [H] shows that 
it is possible to estimate the expectation values of all ob- 
servables of a quantum system by measuring only a single 
"universal" observable on an extended Hilbert space. In 
the Hilbert space of the quantum system, the reduced op- 
erator of this "universal observable" constitutes a min- 
imal informationally complete positive operator-valued 
measurement [15 . Du et al. |16j demonstrated an ex- 
perimental example on this. Allahverdyan et al. |17] 
determine the conditions for making this type of mea- 
surement robust by maximizing the determinant of the 
mapping between the quantum state and the measure- 
ment results. 

The possibility of obtaining the full quantum state in- 
formation from a single measurement appears highly at- 
tractive and may well have practical advantages since it 
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FIG. 1: Schematic representation of different schemes for the 
experimental determination of the state of a quantum system 
S. (a) Quantum state tomography: a series of mutually incom- 
patible measurements are performed, projecting the quantum 
state e.g. along the x, y, and z-axes. (b) The present ap- 
proach: Part of the information in the quantum state is first 
transferred to the assistant A. A single measurement of the 
combined system S -I- A, in direction 9, can then determine 
the complete initial state of S. 



avoids some experimental uncertainties related to mea- 
surement setups for incompatible observables. It re- 
quires, however, the redistribution of the information 
within the extended Hilbert space. This is achieved by 
coupling the system S, whose state is to be determined, 
to an assistant A and letting the combined system evolve 
for a suitable period. The sketch of this measurement 
idea is shown in Fig. [l]b). As we show below, the suc- 
cess of the resulting measurements depends on the form 
of the Hamiltonian as well as on the duration of the evo- 
lution and the choice of the final measurement on the 
combined system. 

In this paper, we study the details of this type of mea- 
surements using a (nuclear) spin 1/2 as the system whose 
quantum state is to be determined, and a different spin 
1/2 as the assistant. We consider in detail what types 
of Hamiltonian can be used to couple the system to the 
assistant, how the information content of the resulting 
state can be maximized, and under what conditions the 
scheme will fail. As an experimental example, we present 
results from a nuclear spin system, using nuclear mag- 
netic resonance (NMR). 



II. COUPLING SYSTEM AND ASSISTANT 

A. Hamiltonian 

We consider two qubits interacting with local magnetic 
fields and coupled through the Heisenberg interaction. 
The system Hamiltonian can then be written as 



H — H;,{Bi, B2) + Hex{Jx, Jy, Jz) 

= BiSl + i?2'S'z + 

jxS^Sj. ~\~ jySySy j^S^S^^ 



where = ^a'^{i' = x,y,z) denotes the local spin opera- 
tor for qubit k. The B^s are the strengths of the external 
magnetic fields (along the z axis) acting on qubit k, and 
the Ji,s are the Heisenberg exchange constants. 

For arbitrary J^, this is often called the anisotropic 
Heisenberg XYZ model. Some special cases are: 

• XXX (or isotropic Heisenberg): = Jy = Jz 

• XXZ ■ Jx ^ Jy Jz 

• XY : J,^0 

• XZ: Jy = 

• Heisenberg-Ising : Jx =^ Jy = 

> and < correspond to the antiferromagnetic 
and ferromagnetic cases, respectively. In many solid- 
state systems, the coupling constants Ji, can be tuned by 
external fields and many proposals for solid-state quan- 
tum information processors rely on their tunability. 

The Hamiltonian of Eq ([T]) splits into three mutually 
commuting parts. 



Hn -t- Hn 



where 



Hzz — JzS^S'^ 

Ho = Bjb{SI - 4') + iiS'+Sl + SlSl) 
H2 - B{Sl + Si) + f 7./(^i^+ + SlSl). 



(2) 



(3) 



B— {B1+B2) /I and J — (Jj;+Jj,)/2 are the average field 
and the coupling constant, and 73 — {Bi — B2) / {B1+B2) 
and 7j — {Jx — Jy) / (Jx + Jy) are anisotropy parameters. 
Sj. — S*^' -I- iSy are the raising and lowering operators. 

With this decomposition, the eigenvalues and eigen- 
vectors can be easily calculated by diagonalizing the sub- 
spaces consisting of Hq and H2 ■ We obtain for the eigen- 
values 



Ai = ^Jz + Vi, 
A2 = ~lJz + V2, 
A3 = -^Jz - V2, 
A4 ~ jJz — f]l, 



(4) 



and for the eigenvectors 



Here 



/ cos(6'i/2) 



V sin(0i/2) 



V / 



( \ 

9i 



COS 



V / 

/-sin(0i/2) 



V cos(6li/2) 



(5) 



r;i = ^B^ + (J7j/2)2 
m = ^{BiB? + {J /If 



(6) 
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and 



III. MEASUREMENT PROCEDURE 



cos -TT = 
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1 ^ , ^r. . = sgn{J)J 

y'2ri2{v2+B'-fB) V 



2j7i 



(7) 
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B. Evolution 

We write the evolution operator as a product of the 
evolutions generated by the three mutually commuting 
terms of Eq. (2): 



Uir) 



where 

Mr) 



-iHr 



U^z{r)U^{r)U2{T) 



(8) 



1+COs(t)2"'') - 



1 — cos(r;2'r) 



^2cos6l2 sin(r?2r)(5i - Sl)+ 
i sin 6*2 sm{ri2T){SlSl + SlS'l] 



(9) 



1+cos(77it) l-cos(T)ir) ^2-j 

i2 COS 6*1 sin(r;ir)(S'] + 5"^)+ 
i sin 6*1 sm{T]iT){SlSl + SlS^_). 

In the following, we will use a different operator basis 
for the diagonal terms: we define the polarization oper- 
ators 1°''^ = |l ± si^\ In terms of these operators, the 
total propagator becomes 

U{t) = aj^l^ + 02/"/^ + ag/f + a4/f 

+d{SXS^_^ + SlSl) + b{SXSl + SlSl), (10) 

where 



(11) 



As we show in the following section, the evolution of 
Eq. ([s]) transfers information between qubits in such a 
way that it becomes possible to measure the complete 
quantum state of one qubit with a single apparatus, as 
proposed by AUahverdyan et al.|17j. 
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A. Principle 

Consider a two-level system S (spin-1) whose state can 
Pii P12 ^ 



be represented by p = with the normaliza- 

\ P21 P22 ) 

tion p\\ -|- P22 — 1- To determine the state p, we can 
measure the vector s = 2Tr{Sp) = {s^, Sy, Sz)'^ where 

Si: = Pl2 + P21, Sy = i{pi2 - P2l) , Sz = Pll " P22 and 

S = {Sxi Sy, Sz)'^ ■ 

To transfer part of the state information to the assis- 
tant A, we couple the two subsystems with the interac- 
tion Hamiltonian H of Eq.|l]). At the time t — 0, the 
composite system S+A is in the state 
Without loss of generality, we assume ^ = ^1 -I- eSz, (0 < 
e < 1). Here, the superscripts S and A refer to the two 
subsystems. 

Under the effect of the coupling Hamiltonian of Eq. 
(1), this state evolves into Qr — IJ{t)qqW{t). On 
this state, we perform repeatedly measure the sim- 
plest possible nondegenerate, factorized observable Cl — 
^a=i '^aPai which determines the complete set {Pa} of 
the joint probabilities, They correspond to the eigenval- 
ues of Qr in the eigenbasis of the observable Cl. Since 
these values were generated from the initial state by a 
one-to-one mapping. Pa = Pkq = Yl^tj ^kq,ijPij, we can 
invert this mapping to calculate the original state p |17j . 

The precision of the back-calculation depends on the 
size of the determinant A = det{Mi^q ij): If |A| is small, 
any (experimental) error in the measurement of Pkq will 
result in a large error in , roughly oc 1/|A|. We there- 
fore seek to maximize |A| and thereby the precision of 
the measurement. This maximization is achieved by a 
suitable choice of the Hamiltonian H, the duration r, 
and the single observable f2. 



B. Symmetry properties of the evolution 

The Hamiltonian H of Eq. ^ consists of three com- 
muting parts: Hzz, Hq, and H2. All of these terms are 
invariant under 7r-rotations around the z-axis. We use 
this property to separate the density operator into two 
parts that transform irreducibly under this symmetry op- 
eration: One part which is invariant with respect to tt^ 
rotations, and the second part, which changes sign. The 
first part includes diagonal terms, zero quantum and dou- 
ble quantum coherence [TH]; the second part includes the 
single quantum coherence terms. The symmetry of the 
Hamiltonian implies that the evolution does not transfer 
information from one subspace to the other. Figure [2] 
illustrates the division of the density operator into these 
subspaces. 

A consequence of this separation of the system into two 
distinct subspaces is that it restricts the possible choice of 
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FIG. 2; Subspaces of the density operator that are invariant 
by the evolution under the Hamiltonian of Eq. 



and its determinant is 

A = 83(Ali2 12-A^11.12)(Alll.ll-^12.22-Xi2,ll7Wii,22)- 

(14) 

Here 5(c) denotes the imaginary part of c. Using Eqs. 
Q, Q, pT] ) and (13 1, we find for its absolute value 



|A| = ^\{l~e^)M-Jzr){[sin{20,)sm\f^,T)]^- 
[sin(26'2) sin2(?72r)]2} + 2e[sin(26li) sin^(r;iT) 
+ sin(26l2)sin^(r;2T)]{[l - 2sin2 6*1 sin^irjir)] 
X sin 6*2 sin(2772T) - [1 - 2 sin^ 6*2 sin^ (»72t)] 
X sin 01 sin(2?7iT)}|. 

(15) 



observables. In particular, if we choose the z-components 
of the two spins to the single observable Ct, then all pos- 
sible combinations fall into the subspace that is invariant 
under tt^ rotations and therefore does not provide infor- 
mation about the other subspace. For this paper, we 
choose the x-components of both spins; another, equiva- 
lent choice would be the y-components. 



C. Transfer matrix 

This evolution process, together with the subsequent 
measurement of the single observable Ct, transfers the 
information from the initial state to the set of mea- 



surement results Pa — 
values of the operators 



Pf^q, which are the expectation 



acting on the state gr 



(-l)'=^f))®(il-(-l)^^(^)) 



Pkg^Tr{PkgUiT)goUHT)}. 

The indices k,q can have the values 1, 2. 

We use the transfer matrix A4 to describe this map: 



IV. OPTIMIZATION: MAXIMIZING A 

The size of the determinant |A| of the transfer map- 
ping determines the quality of measurement. Maximizing 
|A| will minimize the statistical error of the estimation 
during the inversion of Eq. (12 1. We can maximize it by 



an appropriate choice of the initial condition of the assis- 
tant, the parameters of the Hamiltonian that generates 
the evolution U, and the duration of the evolution. 

Figure [3] plots the maximum possible determinant size 
|A|,„Qa: for the Hamiltonian of Eq. ([l| as a function of 
the polarization e of the assistant. Clearly, the quality of 
the measurement should increase with increasing polar- 
ization of the assistant. The dashed line in Fig. [3] shows 
for comparison the maximum possible value for a gen- 
eral exchange interaction, taken from Ref. [17 . At the 
extreme cases of zero and full polarization, the Heisen- 
berg coupling Hamiltonian allows one to reach the max- 
imum possible value, but for intermediate polarizations, 
its maximum value is slightly lower than for the general 
case. 

Let us now focus on two extreme situations e — I (a 
pure state) and e = (a completely disordered state). 
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FIG. 3: (Color Online) Maximal determinant size | A|mai: ver- 
sus the polarization e by the Heisenberg exchange interaction 
([1]), compared to the general case of arbitrary exchange inter- 
action (Eq. (9)) in Ref. [17]). 



(13) 
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A. Assistant in pure state 



|A| reaches its maximum of 1/32 when 



When the assistant A starts in a pure state 1^"^' = 
|l + 51"*^ (corresponding to e = 1), the determinant 
becomes 

|A| = i|[sin(2(?i)sin2(7;iT) +sin(202)sin2(7;2T)] 
X {[1 - 2 sin^ 6i sin^(?7iT)] sin 02 sin(2r72r) 
-[1 - 2sin2 6*2 sin2(ry2T)]sin6li sin(2r;iT)}|. 

(16) 

We can see from this expression that | A | is independent of 
the couphng strength along the z-axis. A Heisenberg 
XY interaction is therefore sufficient for optimizing the 
evolution. We therefore speciahze to this case. Using the 
substitutions 



sin ^ = sin 9k sm{r]kT) 

sinAfe ^ cos(ryfcT)/ cos {k = 1,2), 

we rewrite lAI as 



(17) 



|A| = 



16 



(sin Si cos Ai + sin ^2 cos A2) 



X (cos Si sin ^2 sin A2 — cos ^2 sin ^1 sin Ai) 



(18) 



In terms of these parameters, an optimal solution (|A| 
l/(12-v/3)) is given by the following set of parameters: 



Ai = A2 = A 
sin(2A) = ±1 



— 2 

= ±1. 



(19) 



This parameter set corresponds to the following param- 
eters of the Hamiltonian (1): 



i]kT = TOTT ± i arccos(— Ffc) 
5 = ±771 7(1 



± 



ri)/(i + ri) 

(i+ri)(i-r2) »)2 



IB 

J = 



7J r2(i+ri),,2- 



(i-ri)(i+r2) j)i 
±4?72v/2r2/(l + r2) 



(20) 



sin^ ^, (fc = 1,2). 



%/3 1 
6 ' 2 



% or 



Here, m is an integer and 

(ri,r2) take the pairs of values {\ 

(5 + ~ ^)- Without loss of generahty, we set 

T = 1. The optimal Hamiltonian is then 

H°P^ = 1.1458S']-0.2935S'f-l- 3.38205^5'^ -1.2747S'^S'2 

(21) 



B. Completely disordered assistant. 

When the assistant A is initially in a completely dis- 
ordered state = ^1 (e = 0), we have 



3^|sin(-J,r)||{[(sin(2f?i)sin2(,7iT)]2 
-[(sin(2^^2)sin2(r;2r)]2}|. 



(22) 



sin( J^t) = ±1 ^ J^T = -{2n - 1) (23) 

[n integer) and simultaneously 

sin(26'i)sin^(?7iT) = ±1 and sin(26l2) sin^(r;2T) = 

(24) 



or 



sin(26'i)sin2(77iT) = and sin(26l2) sin^(r72T) = ±1. 

(25) 

Condition (24 1 corresponds to the following parame- 



ters for the Hamiltonian: 



\B\r 
\J.-Jy\ 

IB 



\2m-\\^ 

m 

OR + = OR ?72T = llT 



(26) 



and (25) to 



\J\r 
IB1-B2I 
Bi + B2 



\2m~l\^ 

\J\ (27) 
OR 7j OR T]iT = In, 



where m,l are integers. 

These relationships define six classes of Heisenberg ex- 
change interactions that optimally transfer information 
from the system to the combined system plus assistant. 
The transfer is determined by the product of the Hamil- 
tonian and the evolution time r. Without loss of gener- 
ality, we choose t — j. In these units, some possibilities 
are 

(a) XYX model: = V2{Sl+S^J + 2{SlSl + SlSl) + 

2(1 - 2^/2) S*! as shown in Ref.[T7|; 

(b) XXZ model: H°p' = V2{Sl - S'^,) + 2V2{SlSl + 

y y) '^^z^z'j 

(c) XZ model: H°p' = V2{Sl±S'^,) + 4V2S},Sl+2SlSl 



FAILURE ANALYSIS 



The measurement scheme fails when A = 0. From Eq. 



(151, we see that this occurs when 

sin(26'i) sin2(?7ir) + sin(26'2) sin2(772T) = (28) 
or when 

(1 - e2)sin(-J^T)[sin(26'i)sin2(r;ir) - sin(26'2) sin2(772r)] 
-|-2e{[l - 2sin^ 0i sin^(?7iT)] sin 6*2 sin(2?72r) 
-[1 -2sin2 6*2 sin2(r;2T)]sin6'i sin(277iT)} 0, 

(29) 

independent of the initial state of the assistant A. 

A simple case is sin(20i) = sin(202) = 0, i.e., J = or 
B = or (7b = and 77 = 0). These cases correspond, 
e.g., to 



6 



• a weakly-coupled liquid-state NMR Hamiltonian 

(J, = Jy=0) 

• any Heisenberg interaction without external field 

• an isotropic Heisenberg interaction in the XY plane 
in a uniform external field. 

In all of these cases, the resulting evolution cannot 
generate a state that allows one to measure the complete 
information. 



Another case that fulfills Eq. (|28|) is 
sin(26'i) 

and 



sin(26l2 



7J 

IB 



9 



^ (30) 



|sin(77ir)| = |sin(?72r)| ^ r;iT = \rm: ± t]2t\ 
If, e.g., 771 — r/2, we get the condition 

(7B,7J) = (1,±1) or (-1,1) 

for A = 0. 



(31) 



(32) 



From Eq. (15 1, we can seen that for e = 1 (a pure 



state), |A| does not depend on Jj, while for e — (com- 
pletely disordered state), |A| depends strongly on Jz- For 
e = 0, it is obvious that A = when 



sin( J^r) = => JzT 



(33) 



Hence, the existence of the coupling along the z-axis 
{i.e.jJz ^ 0) is a necessary condition for this measure- 
ment scheme when the assistant A is initially prepared 
in a completely disordered state. In this case, the failure 
condition ( 32 1 can be further modified to 



(7B,7J) = (±1,±1), 



(34) 



which means that when any two among J^, Jy^ Bi, B2 are 
equal to zero, the measurement scheme fails. 



VI. QUANTUM SIMULATION OF THE 
EXCHANGE HAMILTONIAN 

In liquid-state NMR systems, the natural Hamiltonian 
for a system of two spins is 



H 



NMR 



= oJiSl + ti'2'5'f + 27r Ji2>S'_^S'f , 



(35) 



where 0^1.2 represent the Larmor angular frequencies 
of the two qubits (in the rotating frame) and J12 the 
spin-spin coupling constant. This is equivalent to the 
Heisenberg-Ising model. As discussed in the preceding 
section, this Hamiltonian cannot be used to transfer the 
information, since the transfer matrix becomes singular, 
A = 0. Therefore, the key to implement this measure- 
ment scheme in a liquid-state NMR system is to first 
perform a quantum simulation of the Hamiltonian ([T|). 
We briefly discuss two techniques for realizing such an 
evolution. 



A. Short period expansion 

Assuming that we can realize parts of the Hamiltonian 
experimentally, we write the total Hamiltonian as a sum, 

L 
k = l 

In general, the different terms do not commute with each 
other, and it is therefore not sufficient to generate them 
sequentially. However, if the evolution under each term is 
sufficiently short, it is possible to approximate the overall 
evolution in this way. 

Using a symmetrized version of the Trotter formula 

^iA+B)r ^ ^(Ar/2)^iBr)^(Ar/2) ^ q(^3) 

we expand the propagator as 



e 

-iH_ 



iHAt 



-iHi 



2 e 



-iH2 



-iHL 



2 e 



-iHL 



-iHi 



^] + 0(At3), (36) 



which approximates the desired evolution to second or- 
der in At. Keeping At short enough, this allows one to 
efficiently simulate the target Hamiltonian ([T]) by con- 
catenating these evolution periods until the correct total 
evolution is reached. 

Our target Hamiltonian can be decomposed into two 



Hzz and H^y 



non-commuting parts Hz 
JySySy. We thus generate the overall evolution isl as 

U{t) = U"\At) = [Uz{^)U,y{At)Uz{^)r + 0{At'), 

(37) 



where r = mAt is the total duration, and 
;7,(^) = e-(^^+^")^, 

and 



UxyiAt) = e 

represent the evolutions under the partial Hamiltonians. 

Taking as an example the XZ model (case (c) in Sec. 
IV B), it is sufficient to choose the number of evolution 
periods m = 2 for t = 7r/4: the resulting approximate 
evolution 

U^'ir) = [UA^)U.y{l)Uz{^)r 



with 



and 



:c.lc.2 



= e-'^iSl+SD^ ^-^Slsl^ ^^(Sl+Sl)^ (gg) 
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has a fidelity of 0.9958 with the target evolution, where 
the fidelity is defined as 



F([/(T),f7"P(r)) 



The Uzife) operator can be implemented by a free evolu- 
tion period under the internal Hamiltonian if we choose 
uji — ±(^2 — \/2vrJi2 and set the duration to di — jq^j^- 

The Uxy{^) operator can be implemented by four ^ 
pulses (corresponding to e~'^^y^) and a free precession 
period of duration d2 = This evolution period im- 

plements e~^^^^'^^ ; we therefore refocus the chemical 
shift terms by inserting refocusing tt pulses in the middle 
of this period. According to Eq. (38), the second set of 
tt/2 pulses should rotate the spins around the — y axis. 
Here, we choose the -l-y-axis instead to compensate for 
the inversion of the axes system by the 7r-pulses. 

The resulting pulse sequence that generates one seg- 
ment of IfPir) is 

r7riir7ri2 d2 , ,i ,,2 ^2 [TriirTriZ 

di — — — TT TT — — — dl, 

I2\y I2\y 2 ^ '-y^ 2 I2\y l2\y ' 

(39) 

where [9]~ denotes a 9 rotation of qubit k around the i> 
axis. 



B. Exact decomposition 

In some cases, it is possible to achieve the exact trans- 
formation by a suitable decomposition of the evolution, 
using, e.g., 



-iRHWr 



For the propagator ([8|, we can use the decomposition 
U = R - e-'"^^-o'' ■ R\ 



(40) 



where H^nag is the diagonal form of the Hamiltonian. The 
transformation 



/cos^ 



R = 



Vsi 



COS 



sin^\ 



(41) 



which diagonalizes the Hamiltonian ([T]), can be imple- 
mented experimentally (up to an irrelevant overall phase 
factor) by the pulse sequence 



"TT" 


1 


"TT" 


2 _ n 


.2. 


-y 


.2. 








1 r i2 






[tt] 


xM-y 


2 









"TT" 


1 


"TT" 


^ 2 


.2. 




.2. 



2 



"TT" 


1 


"TT" 


1 


"TT" 


2 


"TT" 


.2. 


— X 


.2. 


-y 


.2. 


y 


.2. 



(42) 



with Tl 



2|6 



921 



9i > 62 or 6*1 < 02, and i?^^ by the Hermite time-reversed 
sequence. 

The evolution under the diagonal Hamiltonian 



— — and Lp = X or —x tor 



Ud^. 



-iHdiagT 



-iA2T 



\ 



-iXat 



is realized by the pulse sequence 



J 

(43) 



Tl 
2 



T3 


"TT" 


1 


"TT" 


2 


.2. 


— X 


.2. 



"TT" 


1 


"TT" 


2 


.2. 


— a: 


.2. 


— a; 



where rg = ^^^^^j^r, /3i = 

Aj — A2 + A3 — A 



A1+A2 — A3 — A4 

2 



(44) 
-r and 

plements an evolution under the J-coupling alone, the 
second part implements a composite z— rotation of the 
two qubits by angles Pi and ^2- 

An alternative realization of Udiag is achieved by let- 
ting the system evolve under a constant Hamiltonian 
with 



Ai + A2 — A3 



A4 



and 



L02 



T3 



Ai — A2 + A3 — A4 T 
2 ^3 



for a period t^. 

For the example (c) in Sec. IV B, we choose the pa- 
rameters 



1 



Tl 



T3 



4Ji; 



-,T2 



4Ji 



7r(2 + V2) n{2 - V2) 
Pi = ] ,P2= -. ■ (45) 



When Bi = B2, tp = x in the sequence (42 1; when Bi = 
-B2, if = -X. 



VII. EXPERIMENTAL IMPLEMENTATION 

Experiments were performed at room temperature on 
a Bruker Avance II 500 MHz spectrometer equipped with 
a Triple-Broadband-Observe(TBO) probe at the frequen- 
cies 500.23 MHz for ^H and 125.13 MHz for ^^C. For the 
qubit system, we chose ^'^C-labelled chloroform diluted in 
acetone-dg. The "unknown" state p was prepared on the 
spin of the proton nuclei (^H), which served as the quan- 
tum system S (qubit 1), and the spin of the ^'^C nuclei 
was taken as the assistant A (qubit 2). The spin-spin 
coupling constant is J12 ~ 21A.95Hz. The relaxation 
times were Ti = 16.5 s and T2 = 6.9 s for the proton, 
and Tl = 21.2 s and T2 = 0.35 s for the carbon nuclei. 
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A. Experimental procedure 

There are three steps to implement the measurement 
scheme stated above: (?) Preparation of the initial sys- 
tem, (ii) Quantum simulation and (iii) Measurement. 

Any qubit state p = ^1 + s ■ S can be parameterized 
as a vector in the Bloch sphere: 

= (r sin 9 cos (f>, r sin 9 cos 4>, r cos 9)^ , (46) 



1 for a pure state, and < r < 1 
and (j) are, respectively, the polar 



where the amplitude r 
for mixed states and 9, 
and azimuthal angles. 

The combined system was initialized in the state qq — 
p^^') (K) In our demonstration experiment we chose 

a completely disordered state if — ^l'-"^-' that is exper- 
imentally easy to prepare. Such an initial state go — 
pf"^) ® ll^"^^ prepared by the NMR pulse sequence: 



[arccos (r)]^ 



<i>+-!r/2 



(47) 



The first two RF pulses define the amount of spin po- 
larization on the two qubits. The field gradient pulse 
Gz dephases transverse magnetization to eliminate off- 
diagonal terms in the density operator. The last RF 
pulse turns the remaining (longitudinal) magnetization 
of qubit 1 into the desired orientation. The result of 
the preparation was checked using the standard method 
of state determination based on three noncommutative 
measurements of the system S , i.e., ai^\ a-y^^ and ai^\ 
The experimental results are plotted in Fig. [sjc) and d) 
for r — 1 and |. The experimental average fidelity is 
above 0.99. 

For the coupling Hamiltonian that transfers the infor- 
mation from S to A, we chose H°p* of example (c) (sec- 
tion IV B). We performed two different methods for the 
simulation of this propagator ([s]): the "short period ex- 
pansion" (see section VI A) , using the sequence ( 39 1 with 
m = 2, and the exact decomposition (40) (see section VI 
B) with the NMR pulse sequences (|42|) and ([44l. 

After the coupling evolution, we measured the x com- 
ponents of the two spins to obtain the joint probabilities 
Pkq. For this purpose, we rotated the spins to the z-axis, 

using a [f]^'^ pulse and destroyed off-diagonal elements 
by a magnetic field gradient pulse Gz- The populations 
could then be measured by applying another rf pulse to 
each of the spins and measuring their free induction de- 
cays (FIDs). If the two spins are different isotopes, as 
in our case, their FIDs usually have to be measured in 
separate experiments. 

The resulting pulse sequence for the readout is thus 



'TT' 


1.2 


-Gz- 


'TT' 


i 










.2. 


-V 


.2. 


y 



FID, 



(48) 



where i = 1 or 2 denotes qubit i. The measured FIDs 
along with the normalization condition (X^/cq = 1) al- 



lowed us to reconstruct the four diagonal elements (pop- 
ulations) in the density matrix, which correspond to four 
joint probabilities P^q- The information about the state 
p was then obtained by the inverse mapping . 



B. Experimental Results 



Carbon Signals 



Proton Signals 



(a) 




100 -100 

NMR frequency 



100 -100 

NMR frequency 



(b) 



100 -100 

NMR frequency 



100 -100 

NMR frequency 



(c) 



100 -100 

NMR frequency 



100 -100 

NMR frequency 



FIG. 4: Experimental NMR spectra of carbon and proton for 
different initial conditions. </> = and (a) 9 = Q, (b) 9 = ii/2 
and (c) 9 = n. The y-axis denotes the signal amplitude in 
arbitrary units. 

Figure |4] shows the experimentally observed NMR sig- 
nals after Fourier transformation of the corresponding 
FIDs for the proton and carbon spins for the follow- 
ing initial states: (a) 6* = (i.e. p — hi + Sz), (b) 



0,9 ^ 7r/2 (i.e., p 



Sx), (c) 



TT (i.e.. 



/5= ^1-Sz). 



The amplitudes of the different resonance lines corre- 
spond directly to population differences: 



<S'ArM_R(Proton) 
S'ArAffi(Carbon) 



P 



P: 



Ml 



P 



2m 



fj.2, 



(49) 



where p = I for the resonance line with positive fre- 
quency and p = 2 for the negative frequency line. From 
these populations, we determine the initial condition by 
inverting Eq. (12 1. 

Fig. [5] summarizes these results for a series of similar 
experiments, where we chose initial conditions s (r, 9, (j)) 
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pulse sequence was obtained by combining the sequences 
( 42 1 with ( 44 1 . The results that we obtained were simi- 
lar to those represented in Fig. |5] but the fidelities were 
slightly lower. This difference is probably due to the 
larger number of pulses in this experiment. 



C. Precision of the Measurement 

An alternative measure of the precision of the mea- 
surement is the distance D between the experimentally 
determined state Sexp and the "true" input state s. In 



terms of the parametrization (46), the trace distance be 



tween the two density operators is 
^1 



^exp 



Asl (50) 



where As^ — si^^^-* — s^ih' — x,y,z). 

Writing AP for the experimental errors and using the 
definition ( 12 1 of the transfer matrix, we find for the dis- 
tance 



Dis,xp,s)^ E\AP\, 



(51) 



FIG. 5: (Color Online) Experimental quantum stat e to mog- 
raphy for the general initial state s{r, 6, <j)) [see Eq. (46 1]. We 



compare the results from measuring a single observable of the 
combined system S+A (rows (a) and (b)) with the results 
from the conventional measurement scheme using three non- 
commutative measurements of the system S (rows (c) and 
(d)). In both cases, the expectation values Sx,Sy,Sz, are 
shown from left to right as functions of the angles 6 and (p: 
(a) and (c) for a pure state (r — 1) and (b) and (d) for a 
partially mixed state with r = ^. 



where 



and 



E 



Ex — 

Ey = 

E, = 



El 



El 



El 



As, 
AP 


_ 1 
det(M) 


Et 


.1 Afc2 


As, 
AP 


_ 1 
det(M) 


Et 




As, 
AP 


_ 1 
det(M) 


Et 


,1 



(52) 



(53) 



varying 9 from to vr in increments of 7r/8, and (p from 
to 2ti with an increment of 7r/12. In Fig. we show the 

measured components si"'''''' , s'y^^^ , sf^"^ ior pure states 
(r = 1), while (b) shows the corresponding results for 
mixed states with r = The experiments cover a wide 
range of points on and within the Bloch sphere. The 
experimental results clearly show the expected cosine and 
sine modulations ( 46 ) , indicating that the measurement 



The Aki are the cofactors of the minors Mkj of the trans 



network is effective for all these input states. The average 
fidelity over all iV = 9 x 13 measured states is 



p — 



1 ^ 

-Y 



0.99 



for both cases, r — 1 and r — ^. 

The experimental data shown in Fig. |5]were obtained 
with the "short period expansion" technique of section 
VI A, i.e., the propagator U{t) was approximately real- 
ized by Eq. (37 1 by repeating the sequence (39) twice. 



We also repeated the experiment with the "exact decom- 
position" technique. The propagator U{t) was realized 
by the exact decomposition Eq. ( 40 1 . The corresponding 



fer matrix 




The determinant of this matrix is det{Ai) = —jiA. 
Therefore, the error propagation coefficients Ea depend 
only on the mapping M. The smaller they are, the higher 
the precision of the resulting measurement. 

As Eq. (53) shows, the error propagation scales in- 
versely with the determinant A of the transfer matrix 
A4. We illustrate this dependence in Fig. [6^), where we 
plot the two quantities as a function of the coupling evo- 
lution time r. The minima of E occur near the maxima 
of |A|, and when |A| = 0, E tends to infinity. In this 
range, it is impossible to determine the state p by such 
a measurement. A closer look shows that the minima 
of E do not occur exactly at the maxima of |A|. The 
difference arises from the numerators in ( 53 1 . 



In our experiment, the experimental uncertainties are 
AP « 5%. For the chosen experimental parameters, this 
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results in an average distance Davisexp^s) = 0.04 for 
r = 1 and 0.03 for r = ^. The distance measurement D 
and the fidelity measurement F are related by 1 — < 
D < Vl-F^ EH]. 



(a) 




Evolution timeT 

FIG. 6; (Color Online) (a) Error coefficient E (solid line), 
magnitude of the determinant jA| (dashed line) and their 
product (dotted line) vs. evolution time r under the Hamil- 
tonian H"''* of the example (c) in Sec. IV B. (b) The con- 
currence C (i.e, entanglement) evolves with the evolution 
time r under the same Hamiltonian for the different initial 
states p = |(1 + r sinO cos (j)Sx + r sinO cos (j>Sy + rcosOSz) 
with pure states of r = 1,6 — 7r/2,<^ = (red solid line), 
r = 1,0 ^ n/2, = 7r/2 (blue dashed line), r = 1,6 = 0,(1) = 
(black dotted line), and a specific mixed state of r = 0.8, 6 = 
7r/4, (j) = Tr/6 (red dash-dotted line) 



D. Entanglement 

The evolution that transfers information from the sys- 
tem to the assistant can entangle the two qubits with 
each other. In Fig. |6] we quantify the entanglement gen- 
erated and relate it to the precision of the measurement. 
Fig. |6] b) shows the concurrence C during the coupling 
evolution, calculated as 

C{t) = max{xi - X2 - X3 - X4, 0}, 

where Xi(* — 17 2,3,4) are the square roots of the eigen- 
values of 

in decreasing order, and 

is the instantaneous density operator. 



If the initial state is in the xy-plane, the entangle- 
ment between the system and assistant is maximized at 
roughly the same time as the information transfer for 
these measurements is optimized (as quantified by |A|). 
However, for initial conditions oriented along the z-axis, 
the entanglement generated by the specific Hamiltonian 
shows a relatively complicated time dependence and little 
correlation with the precision of the measurement. For 
evolution times close to r « 57r/4, e.g., the entanglement 
vanishes, while the measurement error is minimized. 

The dash-dotted curve in Fig. ^jp) shows the entangle- 
ment that is generated for a partially mixed input state 
with a general orientation (r = 0.8, 9 = 7r/4, (j) = n/d). In 
this case, the concurrence remains below 0.2 and reaches 
zero even at the times where the measurement precision is 
optimized. If the amplitude r is reduced further, the en- 
tanglement vanishes, C'{gt) = 0, but the precision of the 
measurement is not affected. We conclude that entan- 
glement between system and assistant is not an essential 
criterion for the success of this measurement scheme. 



VIII. CONCLUSION 

We have experimentally demonstrated how the com- 
plete state of a quantum system can be obtained from 
the results of repeated measurements with a single, fac- 
torized observable f2. The procedure, which involves a 
controlled interaction between the system under test and 
a second quantum system, was proposed by AUahverdyan 
et al. [T7]. 

In our experiment, we used a Heisenberg-coupling to 
transfer information from the system to the assistant. 
Interactions of this type are found in many physical sys- 
tems: apart from nuclear spins (like in this work), they 
also occur in quantum dots (211 122] , donor atoms in sili- 
con [231 [21] ; quantum Hall systems [53] and electrons on 
helium . 

The precision of this type of measurements depends 
strongly on the details of the interaction between system 
and assistant, on the type of Hamiltonian as well as on 
the duration of the interaction. This can be understood 
by considering the transfer of information from the state 
of the system to the measurement results from the single 
observable: If we describe this transfer of information 
from n elements of the density operator of the input state 
by a matrix M, the rank of this matrix must be n, i.e. its 
inverse must exist. In practice, it is necessary to choose 
a transfer matrix that is far from the singular case, to 
maximize the precision with which the input state can 
be calculated from the measurement results. 

This initial work has demonstrated the basic possibil- 
ity of implementing such measurements on the simplest 
possible quantum system (a single spin 1/2). Of course it 
is possible to extend the scheme to systems of arbitrary 
size. Work in this direction is currently under way. 
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